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1. INTRODUCTION 

Let us consider the mathematical problem [8-10] 

' A^ + Ap(x)^ = 0, xeO, 

an an 

where A denotes the Laplacian, A is the spectral parameter, O is a bounded and simply 
connected domain in R" for any positive integer n with the sufficiently smooth boundary 
90, n is the outward unit normal to the boundary 90, and the positive quantity p(x) 
is assumed to be 1 outside O. Those A-values for which there are nontrivial solutions \1/ 
and \l/o to (1.1) are known as the transmission eigenvalues for (1.1). Obviously, A = is 
always a transmission eigenvalue, which we view as the trivial one. The problem described 



in (1.1) arises in acoustic scattering from a bounded region with 1/a//o(x) denoting the 
wave speed v and also in electromagnetic scattering from a bounded nonhomogeneity with 



refractive index a//o(x) as a function of location. 

The relevant direct problem involves the determination of the transmission eigenvalues 
when the nonhomogeneity p is known. The relevant inverse problem is to determine the 
nonhomogeneity everywhere in the given domain O by using an appropriate set of A-values 
related to the transmission eigenvalues of (1.1). It is already known that the transmission 
eigenvalues for (1.1) can be determined from some far-field measurements [4,5]. Letting 
A = /c^, we see that, corresponding to each transmission eigenvalue A, we have two /c-values, 
namely k and —k. 

The research field of direct and inverse problems involving transmission eigenvalues is 
now very active, and the related literature is growing rapidly, and hence it is impossible to 
provide a complete bibliography on the general topic of transmission eigenvalues. We refer 
the reader to [10] and the references therein to trace the important developments in the 
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field and to [10,15-17] and the references therein for the related inverse problem of recovery 
of a bounded nonhomogeneity from an appropriate set of transmission eigenvalues. 

We are interested in (1.1) in the special case where O is a sphere of radius b centered 
at the origin in R"^ and p(x) is spherically symmetric, which we write as p{x) with x := |x|. 
We further consider a subset of the transmission eigenvalues for which the corresponding 
eigenfunctions are also spherically symmetric. As in [2] we refer to such A-eigenvalues as 
special transmission eigenvalues of (1.1). 

Let us assume that p belongs to the admissible class A, by which we mean that p{x) 
for X G [0, +oo) is positive, continuously differentiable, and equal to 1 for x > 6, and that 
p"{x) exists almost everywhere, where the prime denotes the x-derivative. 

Let us introduce our key quantity 

D{k) := ^}^^ (f)'(h- k) - cos(A;6) (/)(6; k), (1.2) 



where 4>{x] k) is the unique solution to 

(j)" + /c^ p{x) (/) = 0, < a; < 6, 

(/>(0) = 0, (/>'(0) = 1. 
Note that D{—k) = D(k), and hence D{k) is actually a function of /c^. 



(1.3) 



The following result is known [2], and here we restate it in terms of k rather than A 
for later use. 

Theorem 1.1 Consider the special case of (1.1) with O being the three-dimensional ball of 
radius b centered at the origin, where only spherically- symmetric wavefunctions are allowed 
and it is assumed that such wavefunctions are continuous in the closure ofVt. Suppose that p 
belongs to the admissible class A. Then, the corresponding special transmission eigenvalues 
of (1.1) coincide with the k"^ -values related to the zeros of the quantity D{k) defined in (1.2). 

As in [2], with each nonzero special transmission eigenvalue /c^ of (1.1), we associate 
a multiplicity, which is the same as the multiplicity of kn as a zero of D{k). When p 
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belongs to the admissible class A, it is known [2] that D(k) is entire in /c^ and has the 
representation 

D{k) = 7 E{k), E{k) := e" n (l - ^) , (1.4) 

where 7 is a real constant, /c^ for n G N correspond to the nonzero transmission eigenvalues, 
and d is the multiplicity of the trivial zero transmission eigenvalue. Here we use N := 
{1, 2, 3, . . . }. It is known [2] that some /c^ may be repeated, d is at least 1, and the actual 
value of d is determined by p. In the trivial case p{x) = 1, we have 7 = and hence 
D{k) = Q. 

In analyzing a typical transmission eigenvalue problem, one has to deal with a non- 
selfadjoint eigenvalue problem [2,8-10] and hence in general we cannot expect that all 
eigenvalues will be real. The nonselfadjointness in a natural way forces us to consider 
complex transmission eigenvalues and also transmission eigenvalues with multiplicities. 
The mathematical necessity of including complex and nonsimple transmission eigenvalues 
certainly complicates the analysis. The experimental determination of such transmission 
eigenvalues, especially in the presence of nonsimple eigenvalues, presents a challenge and it 
is an interesting and important open question how to measure them and how to determine 
their multiplicities. 

The constant a defined as 

a := dx \/p{x), (1.5) 

Jo 

has the physical interpretation as the travel time for the wave to move from x = to x = b. 

In [2] we have presented the following uniqueness results regarding the determination of p 

corresponding to the special transmission eigenvalues of (1.1) in the spherically symmetric 

case. Note that it is assumed that we know the value of 6, which is a reasonable assumption 

as far as the applications are concerned. 

Theorem 1.2 Consider the special case of (1.1) with O being the three-dimensional ball 
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of radius b centered at the origin, where only spherically- symmetric wavefunctions are 
allowed and it is assumed that such wavefunctions are continuous in the closure of O. 
Suppose that our data set consists of the corresponding special transmission eigenvalues 
with their multiplicities, and assume that there exists at least one corresponding p in the 
admissible class A. Let a be the constant defined in (1.5). We have the following: 

(a) If a < b, then our data set uniquely determines p; in other words, if both pi and p2 
correspond to our data, then we must have pi = P2- 

(b) If a = b, then our data set along with the value of 7 appearing in (I.4) uniquely 
determines p; in other words, if both pi and p2 correspond to our data and to the 
same 7, then we must have pi = p2- 

In this paper we give an alternate proof of Theorem 1.2, a proof different from that 
given in [2], by providing an algorithm to reconstruct p from the relevant data set. Our pa- 
per is organized as follows. In Section 2 we present some preliminary results that are needed 
later on; the key result in Theorem 2.3 is crucial for the unique reconstruction of p given 
in Sections 3 and 4. In Section 3 we present the alternate proof of Theorem 1.2(a) and the 
reconstruction when a < 6. In Section 4 we present the alternate proof of Theorem 1.2(b) 
and the reconstruction when a = 6. In Section 5 we consider the analogous problem for the 
Schrodinger equation. We present Theorem 5.2, which is the analog of Theorem 2.3 and 
which plays a key role in the reconstruction of the potential in the Schrodinger equation. 
We then give an alternate proof of the uniqueness result of Theorem 5.4, a proof different 
from that given in [2], by providing a reconstruction procedure for the potential in terms of 
the data set consisting of the corresponding transmission eigenvalues with their multiplic- 
ities and the parameter 7 appearing in (5.15). In Section 5 we also provide an illustrative 
example showing that we cannot have uniqueness if the data set does not include 7; let us 
mention, though, that the reconstructed potential is outside the admissible class of poten- 
tials A considered in Section 5 of our paper. Finally, in Section 6 we present some explicit 



examples where we display p{x) and the corresponding quantities 7 and E{k) appearing 
in (1.4); however, in each of those examples, p'{x) has a jump discontinuity and hence is 
outside the admissible class A. In one of the examples presented in Section 6 it is shown 
that the same E{k) yields two distinct p{x) quantities, for one of which we have a < b and 
for the other we have b > a. 

In the recovery algorithms given in Sections 3-5, we solve some basic Riemann-Hilbert 
problems and use basic facts related to their unique solutions. For the benefit of the 
readers who are unfamiliar with the theory of Riemann-Hilbert problems, we summarize 
below their formulation and their unique solutions relevant to our paper. Let us use C for 



the complex plane, C"*" for the open upper-half complex plane, C+ for C+ U R, C for 



the open lower-half complex plane, and C for C U R. 

The idea behind solving a basic Riemann-Hilbert problem is to determine a sectionally 
analytic function on C by determining its sections on C"*" and on C~, respectively, from its 
jump value on the real axis R. Mathematically, we need to solve the functional equation 

F{k)-F{-k) = G{k), ken, (1.6) 

where G{k) is relevant only for real values of k and it indicates the jump. In other words, 
given G{k) for /c G R, we need to determine F{k) for k G C"*" in such a way that F{k) 



is analytic in /c G C"*", continuous in /c G C+, and 0{l/k) as /c — )> 00 in C+. Certainly, 



we then have F[—k) analytic in /c G C , continuous in /c G C , and 0{l/k) as /c — t- 00 



in C~. In general G{k) may not have any extension off the real axis, but even if it does 
only the values of G{k) for /c G R are relevant and needed in solving (1.6). For the unique 
solvability of (1.6) it is sufficient to assume that G{k) behaves like 0{l/k) as /c — ?■ ±00 
on the real axis R and that G{k) is Holder continuous on R with a positive index a. The 
latter condition is expressed as 

\G{ki) - G{k2)\ < C\ki - k2r , /ci,A;2GR, (1.7) 
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for some positive constant C independent of k. In the special case a = 1 the condition 
given in (1.7) is known as the Lipschitz continuity of G{k) on R. In fact, in our paper 
the relevant G{k) is Lipschitz continuous on R due to the fact that G{k) has an analytic 
continuation from /c G R to /c G C; even though the relevant G{k) is bounded on R and 
decays as 0{l/k) as fc — )■ ±00 on R, in general G{k) grows exponentially as /c — )■ 00 in 0+ 
and in C~, but as already stated that unboundedness off the real axis is irrelevant in the 
analysis of the Riemann-Hilbert problem given in (1.6). Under the two aforementioned 
sufficiency conditions on G{k) for /c G R, the Riemann-Hilbert problem given in (1.6) is 



uniquely solvable and the analytic section F{k) defined on /c G C+ is explicitly expressed 
in terms of the values of G{k) known for /c G R as 

F{k) = ^ r dt f^\ keC^, (1.8) 

2m y_oo t — k — i(J^ 

where one can ignore iO"*" in (1.7) if A; G C+ and one can interpret the presence of zO"*" by 
stating that F{k) for /c G R must be obtained by first evaluating the integral in (1.8) when 
k G C"*" and then letting k approach its real value from C+. When G{k) satisfies the two 
relevant conditions on R, it may sometimes be possible to solve (1.6) readily by finding a 



familiar function F{k) satisfying (1.6) with the appropriate properties on C+. Since the 
existence and uniqueness are ensured, we can then conclude that that function F{k) must 
satisfy (1.8). We refer the reader to [11,13] for further information on Riemann-Hilbert 
problems and their solutions. 

In our paper the Riemann-Hilbert problem stated in (1.6) arises in (3.4), (4.5), and 
(5.14). For example, the basic idea behind solving the Riemann-Hilbert problem in (4.5) 
is to split 2ikD{k) into the two pieces Q{k) and Q{—k), where we know that 2ikD{k) has 
the behavior 0{l/k) as k — )> ±00 on R and is Lipschitz continuous for /c G R. In our 
specific case it turns out that 2ikD{k) has an analytic extension in k to the entire complex 
plane C with an exponential growth as /c — > 00 in C except when k — > ±00 on the real 



axis. The splitting is such that Q{k) is analytic in /c G C"*" and bounded in C+ and in fact 
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Q{k) = 0{l/k) as /c — 7- oo in C+. Similarly, the other piece Q{—k) is analytic in /c G C" 



and bounded in C and in fact Q{—k) = 0{l/k) as /c — )■ oo in C . 

In Sections 3 and 4, the cases a < b and a = b, respectively, are analyzed, where a and 
b are the constants appearing in (1.5). In both cases we show that the relevant G{k) in (1.6) 
can be split into the appropriate functions F{k) and F{—k) in such a way that we are able 
to recognize F{k) and express it explicitly in terms of a familiar spectral function, and hence 
we are able to reconstruct the bounded nonhomogeneity from the explicitly constructed 
F{k). On the other hand, when a > b, even though we know that the relevant Riemann- 
Hilbert given in (1.6) is uniquely solvable, we are unable to express the corresponding 
F(k) explicitly in terms of a familiar spectral function yielding the nonhomogeneity. Thus, 
our reconstruction method exploiting the relevant Riemann-Hilbert given in (1.6) does not 
seem to yield the nonhomogeneity when a > 6. In other words, the case a > 6 is an open 
problem, and it is not known in that case whether the nonhomogeneity can be recovered 
by a method similar to that used in Sections 3 and 4. 

2. PRELIMINARIES 

We consider the extension of the differential equation in (1.3) to the half line R"*" := 
(0, +oo), namely 

tjj" + k'^p{x)ij = 0, XGR+, (2.1) 

where p{x) belongs to the admissible class A and hence p{x) = 1 for x > b. Let us define 
the travel-time coordinate y as 

y{x):= ds^/p(s), xG[0,+cx)), (2.2) 

Jo 

and note that (1.5) and (2.2) imply that a = y{b). We remark that (2.2) yields 

y{x) = X + a — b, x > b. 



Let f{x;k) denote the Jost solution to (2.1), i.e. f{x;k) satisfies (2.1) and 

f{x;k) =€''''', f\x;k) = ike''''', x>b. (2.3) 

We also note that, when k — 0, (2.1) reduces to •0"(a;; 0) = 0, and hence with the help of 
(2.3) we get 

/(x;0) = l, xg[0,+oo). (2.4) 

Via a Liouville transformation, (2.1) can be transformed into a Schrodinger equation. 
In other words, if we let 

f{y;k) := [p{x)]'/' e-^'^^'-'^^ /(x;fc), (2.5) 

where y is related to x as in (2.2), then f{y; k) becomes the Jost solution to the Schrodinger 
equation 

r (y; k) + ehy- k) = V{y) f{y- k), ye R+, (2.6) 

where the prime now denotes the y-derivative and we have 

f{y-k) = e^^y, ny-k) = ike^^y, y > a, (2.7) 

With the help of (2.2) and a multiple use of the chain rule in taking the derivatives on 
the right-hand side of (2.8), we write (2.8) as the second-order linear differential equation 



d'' [[p{x{y))f^^ 

' V{y)[p{x{y))Y/\ yGR+. (2.9) 



dy^ 
Note that p{x{y)) satisfies 

d[p{x{y))] 



p(x(y))\ = 1, 



y=a 



and hence 

,1/41 ^1 dipJMM'" 



w-(«))i' i,^« -. ,, 



= 0. 

y=a 



As we shall see, (2.9) will be useful in constructing explicit illustrative examples of p{x) 
and V{y) from some appropriate sets of data. 

From (2.4) and (2.5) we conclude the following result. 

Corollary 2.1 Assume that p{x) belongs to the admissible class A. Let f{y;k) be the 
corresponding Jost solution to (2.6) with Viy) as in (2.8). Then, we have 

f{y; 0) = [p{x{y))]'/\ y G [0, +oo). (2.10) 

Let us remark that, when p(x) is in the admissible class A, we have [6,12,14,18] 

f{x;-k) = f{x;ky, ke-R, (2.11) 

where the asterisk denotes complex conjugation. The result in (2.11) will be useful in es- 
tablishing (2.18) and in the reconstruction of p{x) from a data set containing the associated 
transmission eigenvalues. 

The results in the following proposition are already known but we state them with a 
brief proof for the convenience of the reader. 

Proposition 2.2 Assume that p{x) belongs to the admissible class A. Let V{y) be the 
potential obtained from p{x) as in (2.8) and let f{y:,k) be the corresponding Jost solution 
to (2.6) satisfying (2.7). Then, we have the following: 

(a) The potential V{y) belongs to the admissible class A described in Section 5. Conse- 
quently, the Jost solution f{y;k) has the properties outlined in Proposition 5.1. 

(b) The differential equation (2.1) with the Dirichlet boundary condition V'(O) = 0, cannot 
have, for any negative value of k'^, any solutions that are square integrable in x E R"*". 
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(c) The corresponding half-line Schrodinger equation (2.6) with the Dirichlet boundary 
condition at y = cannot have any nontrivial solutions that are square integrable 
in y E R+. Hence, the corresponding Schrodinger operator has no bound states, and 



therefore /(O; k) cannot vanish for k G C+ \ {0}. 

PROOF: When p is in the admissible class A, the corresponding V{y), as seen from (2.2) 
and (2.8), is real valued, compactly supported, and integrable. Thus, V{y) belongs to 
the admissible class A described in Section 5, and (a) is proved. Note that (c) directly 
follows from (b) because of the Liouville transformation given in (2.5), Proposition 5.1, 
and the fact that the bound states correspond to square-integrable solutions of the relevant 
differential equations. Hence, we only need to prove (b). Because of the selfadjointness 
of the corresponding Schrodinger operator, any existing bound states may occur only at 
negative values of k'^. If t(j{x) were a nontrivial square-integrable solution to (2.1) at some 
negative value of /c^, i.e. if t/j E L^(R+), then from (2.1) we would get t/j" G L^(R"'"), and 
hence, e.g. via Fourier transforms, we would have i/j' G L^(R"'"). By the Cauchy-Schwarz 
inequality, we would then have i/ji/j' G L^(R"'"). However, that would imply the existence 
of a sequence Xn converging to +oo such that 

\t(j'{Xn)'ll^iXn)\ ^ 0, Xn ^ +00. 

From (2.1) through integration we would then get 

i*„/.///_\ , l2 / ^_ „/_^ i„/./_m2 



dx[ip{x)]*ip"{x) + k^ dxp{x)\ijj{x)\^ = 0. (2.12) 

Jo 

Using the Dirichlet condition V'(O) = and an integration by parts on the first integral in 
(2.12), we would obtain 

[i/'(x„)] >'(x„) - [^(0)]*^'(0+) - / dx\ij'{x)\'^ + k^ dxp{x)\ij{x)\^ = 0. (2.13) 

Jo Jo 

Since V'(O) = 0, from (2.1) it follows that 'i(j"{0~^) = and hence '^'(0"'") is finite. Thus, 
letting Xn —J" +00, from (2.13) we would get 

) /•OO 

dx\tlj'{x)f + k^ dx p{x) \ij{x)f = 0, (2.14) 

Jo 
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which is a contradiction because the left-hand side of (2.14) would be strictly negative due 
to the fact that /c^ < 0, 'i/j{x) is assumed to be a nontrivial solution, and p{x) > for 
X e R+. I 

Let (/){x; k) and (/){y; k) denote the solutions to the initial-value problems on the half 

line that are respectively given by 

(p"{x;k) + k'^p{x)(p{x;k) = 0, x G R+, 

(2.15) 
(/>(0;A;)=0, (/>'(0; A;) = 1, 

4>"{y; k) + k^{y- k) = V{y) 4>{y; k), ye R+, 

(2.16) 
4>{0;k)=0, 4>'{0;k) = i, 

where V{y) is related to p{x) as in (2.8) and x and y are related to each other as in (2.2). 
Note that (2.15) and (2.16) are uniquely solvable [7] and that the corresponding solutions 
are entire in /c^. We remark that (2.15) is actually the extension of (1.3) from the interval 
(0,6) to R+ and hence we use (j){x;k) to denote the unique solution to both (1.3) and 
(2.15). 

The result in the following theorem is crucial for the reconstruction of p{x) from the 
data containing transmission eigenvalues, and it will be used in Sections 3 and 4. 

Theorem 2.3 Let p{x) belongs to the admissible class A. Then, the quantity D{k) defined 
in (1.2) is related to the Jost solution f{x;k) to (2.1) as 

D(k) = ^<°^^>:.{<°^-">. k e C. (2.17) 

Zlrv 

For real k -values, we have 

lin[f{0; k)] = kD{k), ken, (2.18) 

where Im denotes the imaginary part. 

PROOF: Let us express the solution (p{x; k) to (2.15) as a linear combination of the linearly 
independent solutions f{x;k) and f{x;—k) to (2.1), where f{x;k) is the Jost solution 
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(2.19) 



satisfying (2.3). We have 

(/){x;k) f{x;k) f{x;-k) ci{k) 

cP'{x;k)\ ~ [r{x;k) f'{x;-k)\ [c^ik) 

where the coefficients ci{k) and C2{k) are independent of x and are yet to be determined. 
With the help of (2.3) and the second hne of (2.15), we evaluate (2.19) at x = 6 and x = 0, 
respectively, and we obtain 



(/>(6; k) 
<P'{b;k) 



ikb 



ikb 



e e 

f{0;k) f{0;-k) 
r{0;k) f'{0;-k) 



ci{k) 

C2{k) 



ci{k) 

C2{k) 



(2.20) 



(2.21) 



From (2.20) and (2.21), by eliminating ci(/c) and C2{k) we get 



(/>(6; k) 
<P'{b;k) 



pikh 

ike^^^ 



p—ikh 

-ike~'^^^ 



f{0;k) f{0;-k) 
/'(0;fc) f'{0;-k) 



(2.22) 



Let [g; h] := gh' — g'h denote the Wronskian of any two functions g and h. It is known 
[7] and can also directly be verified that the Wronskian of any two solutions to (2.1) is 
independent of x. With the help of (2.1) and (2.3) we get 



[f{x;k);f{x;-k)] 



-2z/c, 



and hence 

V(0;fc) /(0;-fc) 

r(0;fc) /'(0;-fc) 
Using (2.23) in (2.22) we obtain 



-1 

20; 



/'(0;-fc) -/(0;-fc) 
-/'(0;fc) /(0;fc) 



-1 

2ik 



,<P'{b;k) 
Writing (1.2) as the matrix product 



ikb 



ikb 



e e 



-f{0;-k) 
f{0;k) 



D{k) 



cos{kb) 



sm{kb) 



c^{b;k) 
<P'{b;k) 



(2.23) 



(2.24) 



(2.25) 
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and using (2.24) in (2.25), after some simplification we obtain (2.17). Finally, using (2.11) 
in (2.17), we get (2.18). | 

3. RECONSTRUCTION OF p{x) WHEN a<b 

In this section we give a proof of Theorem 1.2(a) by providing a reconstruction algo- 
rithm for the unique recovery of p{x) in terms of the data consisting of the corresponding 
special transmission eigenvalues with their multiplicities. Thus, our data set is equivalent 
to the set of zeros (including the multiplicities of those zeros) of the quantity E{k) defined 
in (1.4). Equivalently, the knowledge of our data is equivalent to the knowledge of E(k). 
Using (2.5) in (2.17), with the help of (1.4), we get 

e^fc(fc-a)/(0;A;)-e-^"(''-")/(0;-A;) 

^^^^ = 2ik^nmF~^ • ^'-'^ 

Note that we assume that a < b, where a is the constant defined in (1.5). For the 
reconstruction, we assume that the existence problem is solved, i.e. we assume the existence 
of at least one p in the admissible class A corresponding to our data. The uniqueness 
aspect in the recovery of p{x) follows from the uniqueness in each of the reconstruction 
steps outlined below: 

(a) When p{x) is in the admissible class A, as stated in Theorem 2.2(a), the corresponding 
Jost solution f{y; k) given in (2.5) satisfies the properties listed in Proposition 5.1, 
and in particular (5.3) holds. Thus, using (5.3) in (3.1) we get the large-/c asymptotics 
of E{k) for /c G R as 

Since (b — a) is assumed to be positive, from (3.2) we can determine the value of (b — a) 
and hence a and also the value of 7 [/o(0)]^/^. Note that we do not have the values of 
7 and [p(O)]^/'^ separately, but as we will see this does not create an obstacle for the 
reconstruction of p{x). 
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(b) Letting 

P{k):=mk)-1, (3.3) 

we write (3.1) as 

where we have defined 

ip{k) ■=2i \kE{k)'y[p{0)]^/^-smk{b-a)] . (3.5) 

By the previous step given in (a) above, we know that our data set uniquely determines 
the value of (^(k). Furthermore, using (3.2) in (3.5) we get </?(/?) = 0{l/k) as /c — )■ ±oo. 
When p is in the admissible class A, it is known that E{k) is entire in k'^ and hence 
also in k. The Lipschitz continuity of ip{k) for /c G R follows from the fact that the 
right hand side in (3.5) has an analytic extension to the entire complex plane and that 
ip{k) = 0{l/k) as /c -)■ ±cx). 

(c) Note that (3.4) constitutes a Riemann-Hilbert problem on the complex plane where 
the function f{k) is specified for /c G R and it satisfies the Lipschitz continuity on 
R and behaves as 0{l/k) as /c — ?■ ±oo on R. The goal is to obtain e'^'^^'^~"''P{k) 
and e~^'^^'^~"'' P{—k) in such a way that e'^'^^'^~"''P{k) is analytic in C"*", continuous 



in C+, and 0{l/k) as A; — )■ oo in C+. Since b > a, those properties of e'^^^'^~°'' P{k) 
follow if and only if P{k) satisfies those properties. The unique solvability of (3.4) 
follows from the Lipschitz continuity of ^{k) on R and the fact that ^{k) = 0{l/k) 
as /c — )■ ±oo. When p in the admissible class, we already know from Proposition 5.1 



that /(O; k) is analytic in /c G C+, is continuous in /c G C+, and satisfies (5.3). Thus, 
the function P{k) given in (3.3) helps us to obtain the unique solution to (3.4). As 
we have indicated in Section 1, the unique solution to the Riemann-Hilbert problem 
given in (3.4) is then expressed with the help of (1.8) as 

,..(.-»,p„)._L£,,_|W_, ,ecF, ,3,6) 
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We see from (3.3) and (3.6) that 



e-'f^i^-^) .00 ^^^^ 



mk)^i^^— it^-f^, fc.c+. 



00 



(d) Having constructed f{0;k) from E{k), we remark that, as stated in Theorem 2.2(c) 
/(O; k) does not have any zeros on the positive imaginary axis in the complex k- 
plane and hence the corresponding half-hne Schrodinger equation with the Dirichlet 
boundary condition has no bound states. Thus, we can use the Marchenko procedure 
described in Proposition 5.1 to uniquely reconstruct V{y) of (2.8) and the Jost solution 
f{y; k) from the already reconstructed /(O; k). This is done by first constructing the 
scattering matrix S{k) from /(O; k) as in (5.5). Next, the Marchenko kernel M(^) is 
constructed as in (5.9) but without the summation term there due to the fact that 
there are no bound states. Then, the Marchenko integral equation (5.8) is uniquely 
solved for K{y,^), and the quantities V{y) and f{y;k) are uniquely recovered as in 
(5.10) and (5.11), respectively. 

(e) Having recovered f{y;k)j we obtain p{x{y)) from f{y;0) with the help of (2.10), 
namely 

p{x{y)) = [f{y;0)]\ ye[0,+oo). (3.7) 

Our next task is to obtain p in terms of x by establishing the relationship between x 
and y. From (2.2) we have 

^ = v%^). (3.8) 

Hence, from (2.2) and (3.8) we obtain the first-order, separable ordinary differential 
equation 

^ dx, xe R+, (3.9) 



with the initial condition y(0) = 0. By integrating (3.9), the relationship between x 
and y is obtained as 

^= / TTT^^' ye[0,+oo). (3.10) 
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Since p{x) is assumed to be positive, from (3.10) it follows that the mapping x \-^ y 
is one-to-one and onto on R"'". Having x as a function of y in (3.10), we can invert it 
to get y as a function of x. Thus, by using (3.10) in (3.7) we recover p{x) in terms of 
X as 

pix) = [fiy{x),0)]\ a;e[0,+oo). 

Thus, the reconstruction of p{x) for x G [0, +oo) from E{k) for /c G R is accomplished. 

Finally, let us note that our procedure yields the value of the constant 7 from the 
knowledge of E{k). This is because we already have the value of 7 [p(O)]^/'^ from (3.2) and 
we have the value of p(0) from (3.7) evaluated at y = 0. 

4. RECONSTRUCTION OF p{x) WHEN a = b 

In this section we consider the case a = b, where a and b are the quantities appearing 
in (1.5). We give an independent proof of Theorem 1.2(b) by providing a reconstruction 
algorithm for p{x) from the data consisting of the corresponding special transmission eigen- 
values with their multiplicities and the constant 7 appearing in (1.4). By Theorem 1.1 the 
knowledge of our data is equivalent to knowing the zeros (with multiplicities) of the quan- 
tity D{k) given in (1.4) as well as the value of 7 there. Hence, our data set is equivalent 
to the knowledge of D{k). 

As seen from (1.4) and (3.2), if D{k) = Oil/k"^) as k ^ ±00, we can deduce that 
a = b. 

When a = b, let us outline the unique recovery of p{x) from D{k) given in (1.4). 

(a) Since a — b, from (1.4) and (3.1) we see that 

2..i,w ^ Z<Ml_il^, ,,c. ,4.:) 

On the other hand, from (2.10) we know that 

/(0;0) = [p(0)]l/^ (4.2) 
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and hence we rewrite (4.1) as 

2«B(,) . /(ft'')-/(ft-^), , , C. (4.3) 

/(0;0) 

We know from (4.2) that /(O; 0) is real and in fact positive because p(0) > 0. When 
p{x) is in the admissible class A, by using (5.3) and (5.4) in (4.1) we conclude that 
2ikD{k) — 0{l/k) as /c — ?> ±oo on R. Furthermore, we know that D(k) has an analytic 
extension to the entire complex plane C and hence 2ikD{k) is Lipschitz continuous 
on R. As seen from (5.3) and (5.4), 2ikD{k) is unbounded as /c — ?■ oo in C"*" and in 
C~, but for the analysis of the Riemann-Hilbert problem to be studied we need the 
laige-k asymptotics of 2ikD{k) only on the real axis. 



(b) Letting 



e(fc) .^ im^. (4,4) 

/(0;0) 



we write (4.3) for real /c-values as 

Q{k)-Q{-k) = 2ikD{k), ken. (4.5) 

Note that (4.5) constitutes a Riemann-Hilbert problem on the complex plane where 
the function 2ikD{k) is specified for /c G R and it is Lipschitz continuous on R and 
behaves as 0{l/k) as /c — )■ ±oo on R. The goal is to obtain Q{k) and Q{—k) in such 



a way that Q{k) is analytic in C"*", continuous in C+, and 0{l/k) as A; — ?> oo in C+. 
Since the unique solvability of (4.5) is assured by the two relevant properties of 2ikD{k) 
stated in (a), we know that the function given in (4.4) must be that unique solution. 
The corresponding properties of Q{k) are deduced from (4.4) by using the relevant 
properties of f{0;k) given in Proposition 5.1. In particular, using Proposition 5.1(a) 
we establish the analyticity of Q{k) in C"*"; using in (4.4) Proposition 5.1(a) and the 



fact that /(O; 0) > we conclude the continuity of Q{k) in C+; using (5.3) in (4.4) we 



obtain Q{k) = 0{l/k) as A; — )■ oo in C+. Therefore, as indicated in (1.8), the unique 

18 



solution to the Riemann-Hilbert problem given in (4.5) satisfies 



oo 



We see from (4.4) and (4.6) that 



/(0;0) /•" ,„ tD{t) 



;(0;fc) = l + ^A_^/ *^-^^. *ec.. (4,7) 



OO 



(c) Let us remark that (4.7) also follows from the Schwarz integral formula for the half 
plane [1] by using the following argument. As stated earlier, Q{k) is analytic in C"*", 



continuous in C+, and 0{l/k) as /c — ?> oo in C+. For real /c-values, from (4.5) and 
(5.2), we obtain 

kD{k)^lm[Q{k)], keK. (4.8) 



Thus, we can construct Q{k) for k E C+ from its imaginary part known for /c G R by 
using the Schwarz integral formula [1] 

Hence, (4.8) and (4.9) yield (4.6) and in turn (4.7). 

(d) Having constructed Q{k) from D{k), we can evaluate the value of /(0;0) by setting 
fc = in (4.4), which yields 

^^<°^°' = d?(o)- <"»' 

Thus, (4.6), (4.7), and (4.10) imply that D{k) uniquely determines f{0;k). 

(e) Having reconstructed /(O; k) from D{k), we can then repeat the reconstruction steps 
(d) and (e) of Section 3 in order to uniquely reconstruct p{x). Thus, the reconstruction 
of p{x) for X G [0, h] with h = a from the data set D{k) is accomplished. 

As stated in [2] , when 6 = a it is an open problem whether the value of the constant 7 
appearing in (1.4) is needed or whether 7 can be determined from the data set consisting 
of E{k) given in (1.4) alone. 
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5. RECONSTRUCTION OF THE POTENTIAL IN THE SCHRODINGER 
EQUATION 

The transmission eigenvalue problem for (1.1) has an analog for the Schrodinger equa- 
tion. The transmission eigenvalues in that case correspond to those values of k'^ for which 
there exist a nontrivial solution pair \1/ and \l/o to the system 

-A^o = k^^Q, y e O, 



(5.1) 



an an 



where V{y) is a real-valued potential that is square integrable on O, and it is assumed 
that V^(y) = outside O. In the spherically-symmetric case in R"^, using V{y) instead of 
V{y) with y := |y|, we define the special transmission eigenvalues of (5.1) as those trans- 
mission eigenvalues for which the corresponding wavefunctions are spherically symmetric 
in addition to V being spherically symmetric. 

We remark that the potential V{y) we use in this section does not necessarily come 
from any function p appearing in (1.1) or (2.1) via the Liouville transformation in (2.8). 
The only assumption we make on V{y) is that it is real valued, compactly supported 
within the interval y G [0,a], and integrable on (0,a). We will say that V{y) belongs to 
the admissible class A if V{y) satisfies those conditions. 

The results given in the following proposition are either known or can easily be proved 
by using the available results [3,6,12,14,18] for the half-line Schrodinger equation by ex- 
ploiting the compact-support property of the potential. We provide a brief proof for the 
convenience of the reader. 

Proposition 5.1 Assume that the potential V{y) belongs to the admissible class A with 
support within the interval [0,a]. We then have the following: 

(a) The corresponding Jost solution f{y; k) has an analytic extension from k E H to the 
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entire complex plane C for each fixed y. Similar to (2.11), we have 

J{y--k) = f{y-k)\ A;gR. (5.2) 

(b) The quantity f{0;k) is nonzero in C+ except perhaps at a finite number of points on 
the positive imaginary axis, say at k = ijSj for j = 1,...,A^ for some nonnegative 
integer N. Such zeros are all simple and they correspond to the bound states of the 
half-line Schrodinger equation with the Dirichlet boundary condition at the origin. 

(c) The large-k asymptotics of f{0;k) in C are obtained via 

f{0-k) = l + o(^j , /c^oomC+. (5.3) 

f(0;-k) = l + o(^^+e-^'''''o(^\ A;^cx)inC+. (5.4) 

Hence, e^'^''/(0; —k) is bounded in C+ for any b satisfying b > a. 

(d) The scattering matrix S{k) for the half-line Schrodinger equation with the Dirichlet 
boundary condition is defined as [3,6,12,14,18] 

and it has a meromorphic extension from k E H to k E C"*" with simple poles occurring 
at k = ijSj for j = 1, . . . ,N. 

(e) The scattering matrix satisfies 

^(/c) = l + o(^J +e-2*'="oQj , A;^ooinC+. (5.6) 

(f) Associated with each bound state there is a positive number, known as the correspond- 
ing norming constant, defined as 

^ j^ dy[f{y,im^ 
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Because V{y) has support confined to the finite interval y G [0,a], the norming con- 
stants are uniquely determined by /(O; k) alone, or equivalently by the scattering matrix 
S{k) alone, as 



c, = y'iRes(^(A;),z/3,), j = l,...,N, (5.7) 

where Kes{S{k),i(3j) denotes the residue of S{k) at the pole k = i(3j. 

(g) The potential V{y) and the Jost solution f{y;k) are reconstructed from the solution 
K{y^ ^) io the Marchenko integral equation 

/•OO 

K{y,0 + M{y + 0+ ds K{y,s) M{s + 0, < y < ^ < +00, (5.8) 



where 

N 



1 r'^ 
M(0 := — / ^^ [1 - S{k)] e'^^ + Y^ c]e-^^^, ^ G R. 



(5.9) 



In fact, we have 



V(v) = -2^?^, V > 0, (6.10) 

/•OO 

f{y;k)=e''y+ dsK{y,s)e''\ (5.11) 

PROOF: The results are mainly known [3,6,12,14,18]; for example, (5.3) and (5.4) can be 
derived by using the integral representation [3,6,12,14,18] for the Jost solution, namely by 
using 

/(y; k) = e^'^y + lfds [sin k{s - y)] V{s) f{s; k). 

We then get (5.6) by using (5.3) and (5.4) in (5.5). The proof of (5.7) can be outlined as 
follows. From (2.7) and (5.11), it follows that K(y,^) = 0f or a<y<^< +00, and hence 
(5.8) in turn implies that M(y + ^) = Ofora<y<^< +00. This fact, combined with (d) 
and (e) allows us to evaluate M{y + ^) for a < y < ^ < +00 by using (5.9) as a contour 
integral along a semicircle in C"*" with its center at the origin and with its radius becoming 
infinite in the limit. Using M{y + ^) = for a < y < ^ < +00 in (5.9), we get (5.7). | 
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Analogously to (1.2), let us define 

D(k) := !l!liM ^'(a; k) - cos(ka) ^(a; k), (5.12) 

k 

where a is the positive constant related to the support of V{y) and (/){y; k) is the unique 
solution to (2.16). The following fundamental result is the analog of Theorem 2.3, and its 
proof is omitted because it is similar to the proof of Theorem 2.3. 

Theorem 5.2 Assume that the potential V{y) belongs to the admissible class A. Then, the 
quantity D{k) defined in (5.12), is related to the Jost solution f{y;k) appearing in (2.6) 
and (2.7) as 

m = m^-m-^), , , c. (5.13) 

ZZK 

For real k -values, we then have 

lm[f{0;k)] = kD{k), keK. (5.14) 

We have [2] the following analog of Theorem 1.1. 

Theorem 5.3 Consider the special case of (5.1) with O being the three-dimensional ball of 
radius a centered at the origin, where only spherically- symmetric wavefunctions are allowed 
and it is assumed that such wavefunctions are continuous in the closure of O. Then, the 
corresponding special transmission eigenvalues of (5.1) coincide with the k^ -values related 
to the zeros of the quantity D{k) defined in (5.12), where (/){y; k) is the unique solution to 
(2.16) with the potential V{y) belonging to the admissible class A. 

When V(y) belongs to the admissible class A, the quantity D{k) defined in (5.12) is 
known [2] to be entire in /c^ and has a representation analogous to (1.4), namely 

D{k)=^e^\{(l-^\ (5.15) 

n=l TT- ' 

with /c^ for n G N being the nonzero transmission eigenvalues, some of which may be 
repeated, and d denoting the multiplicity of the zero transmission eigenvalue. As in [2], 
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we refer to the multiplicity of a nonzero zero k^ of D{k) as the multiplicity of the special 
transmission eigenvalue k'^. 

The following uniqueness result was proved in [2] and is the analog of Theorem 1.2. 

Theorem 5.4 Assume that V{y) belongs to the admissible class A. Then, V{y) is uniquely 
determined by the function D{k) appearing in (5.12) and (5.15) if we assume that there 
exists at least one V{y) in A corresponding to D{k). Equivalently stated, if the existence is 
ensured, V{y) is uniquely determined by the knowledge of the special transmission eigen- 
values of (5.1) with their multiplicities and the constant 7 appearing in (5.15). 

Our goal in this section is to give an independent proof of Theorem 5.4 and further 
provide a reconstruction of V{y) from D{k). The reconstruction consists of the following 
steps and the uniqueness follows as a result of the uniqueness in each reconstruction step. 

(a) First, reconstruct /(O; /c) from D{k), where f{y;k) is the Jost solution appearing in 
(2.6) and (2.7). This is done by solving the Riemann-Hilbert problem given by 

[f{0;k)-l]-[f{0;-k)-l]=2ikD{k), keK, (5.16) 

which is obtained from (5.13). It follows from (5.3) and (5.16) that 2ikD{k) behaves 
as 0{l/k) as /c — )■ ±00 on R. Furthermore, because D{k) has an analytic extension to 
the entire complex plane, it follows that 2ikD{k) satisfies the Lipschitz continuity in 
R. Thus, the Riemann-Hilbert problem in (5.16) has a unique solution that is given 

by 



We can write (5.14) as 

Im[/(0; A;) - 1] = kD{k), keK. (5.18) 

The result given in (5.17) also follows by using the Schwarz integral formula (4.9), by 
replacing Q{k) there with /(O; k) — 1, with the help of (5.18), we obtain (5.17). 
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(b) Having obtained /(O; k) from D{k), we can use (5.5) and (5.7) in (5.9) and obtain the 
Marchenko kernel M(^) from f{0]k). 

(c) The potential V{y) is then uniquely reconstructed as in (5.10) by using M(^) as input 
to the Marchenko integral equation (5.8) and by obtaining K(y, ^) as the unique 
solution to the Marchenko equation. 

Let us mention that it is an open problem whether the value of 7 appearing in (5.3) 
can be determined from the zeros of D{k). If the answer is yes, then 7 is not needed for 
the unique determination of V{y), and the zeros of D{k) with their multiplicities would 
be sufficient for the reconstruction of V{y). In the following example, we show that 7 is 
needed to construct a potential, which is, however, outside the admissible class A. 

Example 5.5 Let the potential V{y) be given as 

V{y) = c6{y-a), 

where c is a real nonzero constant, 6{y — a) denotes the Dirac delta function with argument 
y — a, and a is the positive number related to the interval [0, a] containing the support of 
V{y). The corresponding Jost solution is obtained by solving (2.6) and (2.7), and we get 

hy:k) = \ ^ 2^/cy 2^k (519) 

I e'^'y, y>a. 

From (5.19) we evaluate /(O; k) and then using (5.13) and (5.15), we obtain the values of 
7 and D{k), yielding 

,..„^ ,,„..M.(f!lgf^)^^(,_^)^ ,,.0, 

\ / n=l ^ ^ 

Hence, in this example the transmission eigenvalues, i.e. the /c^-values corresponding to 
the zeros of D{k), all have double multiplicities and are given by /c^ = n^yr^/a^ for n G N. 
However, as seen from (5.20) E{k) alone does not uniquely determine c, and hence c or 
equivalently 7 is also needed for the unique determination of V. 

25 



6. EXAMPLES 

In this section we illustrate the transmission eigenvalue problem corresponding to 
p{x) appearing in (1.3) with some explicit examples. In our first example, with the help of 
(2.2), (2.5), (2.9), and Example 5.5 we present a concrete p{x) for which we can explicitly 
evaluate the relevant quantities D{k) and E{k), given in (1.2) and (1.4), respectively. 

Example 6.1 Let e be a positive parameter and let c be a real nonzero parameter. Assume 
(2.2) is given by 

y{x) 





' 


e^x 


Xo, 






< tcx + 1 






u 


-a:o +?/o, 


x> 


xo. 


e 


-1 


2/0 := y{xQ) -- 


e 


-1 




ec 




c 



where 

Xq:^- -, yQ:=y{xo) = - -. (6.1) 

ec c 

In order to have xq and yo positive, we must have c > if e > 1 and we must have c < 

if e < 1. Note that (6.1) implies that 

(e-l)2 

b- a = xo~yo = , (6.2) 

ec 

where a and b are the parameters appearing in (1.5). We have b > xq and hence from (6.2) 

we see that a < 6 if e > 1 and that a > 6 if e < 1. Using (2.2) we get 

2 f '' 

p(x)=(|) = (ecx + ir ^-^°' (6.3) 

[ 1, X> Xo, 

Because of (6.1), we see from (6.3) that p{x) is continuous at xq, whereas p'{x) jumps from 
p'{xq) = —4c to p'{xq) = 0. One can directly verify that the Jost solution f{x; k) to (2.1) 

is given by 

f ^^e-^'=(-^)VMz(x;fc,e,c), x < x,, 

y e'^^^, X > Xq, 

where we have defined 

Z{x; k, e, c) := {2k + ic)e^'=^'^/^^^^+^) - ice'^''''^'-^^/^-'''''''/^'^''+^\ 

26 



One can check that f{x; k) and f'{x; k) are continuous at x = xq. From (6.4) we get 

1 



/(0;fc) = -^e-^'=(^-i)^/(-) 
2ek 



{2k + ic) —ice 



_ ■ „2ifc(e-l)/c 



Using (6.5) in (2.17) we obtain 



D{k) 



2eA;2 



cos 



k{e-l) 



ec 



cos 



k{e^ - 1) 



ec 



2k . [kit -I) 

— sm 

c \ ec 



By expanding (6.6) in powers of /c^ we get 



D{k) = ^e + o{k^ 



where a comparison with (1.4) reveals that 



d = l, 7 



Using (6.6) and (6.7) in (1.4) we have 



A; ^ in C. 



3e3c3 ■ 



E{k) 



3e^c 



2^4 



2(e-l)4A; 



4^2 



cos 



Me-1) 



ec 



cos 



k{e^ - 1] 



ec 



2/c . fkit-1 

— sm 

c V ec 



From (6.8), we get 
E{k) 



3e2c3 /A;(e-1)2\ /I 



/c — )> ±oo, 



and hence a comparison of (6.9) with (3.2) reveals that 



b — a 



i^-iy 



ec 



7 [p(0)] 



1/4 



3e2c3 



(6.5) 



(6.6) 



(6.7) 



i2\ n 



(6.8) 



(6.9) 



which is compatible with the value of (6 — a) given in (6.2), p{0) from (6.3), and 7 in 
(6.7). Note that p{x) in this example is outside the admissible class A because of the jump 
discontinuity of p'{x) at xq. Let us remark that D{k) given in (6.6) can also be obtained 
by using (1.2), where the unique solution (j){x; k) to (2.15) in this case is given by 



(j){x; k) 



(ecx + 1) . / e'^kx 
sm 



X < Xq, 



e^/c Kecx + 1 

C3{k, e, c) sin(A;x) + C4(/c, e, c) cos(A;a:), x > Xq, 
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with the constants cs{k, e, c) and C4{k, e, c) specified as 



1 //c(l-e)2\ c 

C3(/c,e,c):= — cos^^ +^ 



C4(/i;, e, c) := — - sm i — --— - 



sin 



cos 



k{l-e) 



ec 



k{l- 



€C 



sm 



cos 



A;(l-e2 



ec 



fc(l-e2 



ec 



Using the result of Example 6.1, in the next example we will produce two distinct 
profiles p{x) corresponding to the same E{k) but to different 7 values; in fact, in one case 
we will have a > b and in the other case we will have a < b. 

Example 6.2 In Example 6.1 above, let us use the following values for the parameters 

1 

e = 2, c=-, 



where b is the constant that appears in (1.2) and is related to the known support of p{x) — l. 
Using (6.2), (6.3), (6.7), and (6.8), we obtain 

_16b^ 

p{x) 



7 



(2x + 6)4' 

1, x> 6/2, 

-b^ 36 

= , a = — , 

24 ' 2 ' 



X < 6/2, 



E{k) 



-6 

64p 



bk\ fdbk 

cos I — COS 

2 V 2 



26A; sin ( y") 



On the other hand, in Example 6.1 if we use the parameters 



'=r ' 



from (6.2), (6.3), (6.7), and (6.8), then we obtain 

_64 
p{x 



1 



(26-x)4' 
1, X > b. 



X < b. 



(6.10) 

(6.11) 
(6.12) 



(6.13) 



7 = — 1 a = -, 
'6 2' 



(6.14) 



28 



-6 

E{k) 



, bk\ f?,hk\ ^,, . fhk 

cos — — COS — 2bk sin — 

2 \ 2 V 2 



,.p ...... ^^.. . , ..,....., ^, . (6.15) 

Thus, as seen from (6.12) and (6.15), we have produced two distinct profiles for p{x) given 
in (6.10) and (6.13), respectively, corresponding to the same E{k), but two different 7 
values. In fact, as seen from (6.11) and (6.14), the former corresponds to the case a > b 
and the latter to a < 6. We can simplify and rewrite (6.12) as 

^,,, 12 rbk\ ( sm(bk)\ 

As seen from (6.16), corresponding to the two distinct profiles given in (6.10) and (6.13), 
we have a simple zero transmission eigenvalue, infinitely many simple nonzero real trans- 
mission eigenvalues that are given by /c^ = 4n^7r^/6^ for n G N, and infinitely many simple 
complex transmission eigenvalues that are related to nonzero zeros of kb — sm{kb). Note 
that for each complex transmission eigenvalue, its complex conjugate is also a transmission 
eigenvalue. 

We conclude with another explicit example. 

Example 6.3 For a positive parameter c, let 



{b + c 



i2 



^2' 



X <b. 



p{x)= { {x + cf 

1, X >b, 

where b is the positive parameter appearing in (1.5). Using (2.2) we obtain 

f {b + c) logfl + -) , x<b, 

\ X — b -{- a, X >b, 

where a is the parameter appearing in (1.5) and its value is obtained from (6.17) as 

a:=y{b) = {b + c)\og(l + -\. (6.18) 

In this case a > b because from (6.18) it follows that 
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based on the positivity assumption on c. We can solve (1.3) explicitly and get 



(/)(x; k) 



r+ — r. 



[-c"+(x + c)"-+c"-(x + c)"+], xe [0,6], 



where we have defined 



r± 



1± Vl-4(6 + c)2/c2 



(6.19) 



(6.20) 



Using (6.19) in (1.2) we obtain 



D{k) = ^^4^ <^'(^; k) - cos{kb) (P{b; k), 



(6.21) 



where as we see from (6.19) and (6.20), with the help of r_|_ + r_ = 1, we have 

(f>{h-k) 



r-i- — r_ 



-r^(-r^^ 



ct>'{h;k) 



,1 + - +T-^- 1 + - 
+ c \ c I o + c\ c 



Letting /c — )> in (6.21), with the help of (1.4), we get 



r+ 



rf = 1, 7 = c^ 



2 /6 



3 Vc 



= '')-<-D'-(-^ 



(6.22) 



where, by using a graphical argument, it can be shown that 7 < 0. We remark that the 
results in this example are also valid if c < but h > —c. In that case from (6.18) we get 
a < b and from (6.22) we get 7 > 0. 
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